close all
%clear all
%load Results_baselinecounterf.mat
C3=resultsC3(1,:);
P3=resultsP3(1,:);
gt1=resultsgc3(1,:);
gt3=resultsgc3(1,:);
Chat=(exp(C3));

for i=1:11
    CC1(i)=log(Chat(1))+(i-11)*(gt1(1));
    CC2(i)=log(Chat(TT))+(i-11)*(gt1(150));
end

 for i=12:TT
     CC1(i)=log(Chat(1))+(i+1-12)*(gt1(1));
     CC2(i)=log(Chat(TT))+(i+1-12)*(gt1(150));
 end

CC=log(Chat(1:TT))+cumsum(gt3(1:TT));

test=([(CC1(1:11)) log(Chat(1:TT))+cumsum(gt3(1:TT))])';








%%%%HERE
figure
hold on 
graph_baseline=plot(test(1:100)-CC1(1:100)')
set(graph_baseline,'LineWidth',2, 'LineStyle', '-','color', 'b');
graph_zeros=plot(zeros(100,1))
set(graph_zeros,'LineWidth',2, 'LineStyle', '--','color', 'k');

legend('Baseline','Only IPR', 'Only trade lib.','Initial BGP')
%legend('Baseline','Only IPR','Initial BGP')

%title(['\fontsize{16} China'])
ylabel(['\fontsize{14} log consumption (rel. initial BGP trend)'])
xticks([1 6 11 16 21 26 31 36 41 46 51 56 61 66 71 76 81 86 91 96 101 106 111])
xticklabels({'-10','-5','0','5','10','15','20','25','30','35','40','45','50','55', '70','75','80','85','90','95','100','105','110'})

hold off

